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Abstract 

We study a Heisenberg S=l/2 ring-exchange antiferromagnet which exhibits a quantum phase transition from a spontaneously 
dimerized (valence bond sohd) phase to a magnetically ordered (Neel) phase. We argue that the quantum transition is of uncon- 
ventional nature; both singlet and triplet modes of high density condense as the transition is approached from the dimer side, 
signaling restoration of lattice symmetry. These features are consistent with "deconfined quantum criticality" , of which the present 
model is believed to be the only example so far. 
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Recent work on quantum critical phenomena has re- 
vealed the fascinating possibility of "deconfined critical- 
ity" , where spinon deconfinement occurs at a quantum crit- 
ical point (QCP) separating a valence bond solid (VBS) 
and magnetically ordered (Neel) phase [T] . In order for this 
exotic scenario to take place, the magnetic and VBS or- 
ders must not coexist; otherwise the spinons remain per- 
manently "paired" in S=l quasiparticles |H2j . The intense 
search for a microscopic Heisenberg model that exhibits 
this novel criticality has led to a solid possibility: a S=l/2 
model with four-spin ring exchanges [5] . The Hamiltonian 
of this two-dimensional, square-lattice model reads 



(1) 



i 7 j,k,l 



where the couplings are antiferromagnetic, J, K > 0. The 
four-spin term acts on a given plaquette as — if {(Si -82) (S3 ■ 
S4) + (S2 • Ss)(Si • S4)}, where the spins are numbered as 
in Fig. [2(b). Summation is performed over all plaquettes, 
so that H does not break the square lattices symmetries. 
The quantum Monte Carlo study of Ref. .3] found a critical 
coupling (K/J) c w 2, separating a VBS phase (K > K c ) 
and a Neel phase (K < K c ). The spontaneous dimerization 
in the VBS state was suggested to be of the columnar dimer 
type, as in Fig.QJb). 

In the present work we study the model ([I]) by approach- 
ing the QCP from the dimerized phase. In this case it is 
natural to represent the spins in terms of "bond triplons" 
t ia ,a = x,y,z (3] which create S=l triplet excitations on 
a given dimer i (see Fig. Gib)). These quasiparticles are 



hard-core, meaning that two triplons cannot be created on 
the same site (dimer) . The analysis of the effective Hamil- 
tonian H e a[t a ], written in terms of the triplon operators, 
leads to the spectrum oj(k) which determines the location 
of the QCP where the magnetically ordered state emerges. 

Even though the representation of the original spin 
Hamiltonian in terms of H e ^[t a ] is exact, in practice one 
has to deal with the various interactions between the 
triplons in an approximate way. These interactions in- 
clude both the kinematic hard-core constraint and the 
various dynamic triplon-triplon scattering vertices. For 
the model Jl} up to 8-point vertices appear. The success 
of perturbation theory in treating such interactions is de- 
termined by the quasiparticle density n which it typically 
low (microscopically the density is related to the quantum 
fluctuations). In quantum spin models involving explicit 
dimerization, i.e. situations where some of the exchange 
coupling are stronger than the others (such as spin lad- 
ders, coupled ladders, bi-layer models, etc.), calculations 
based on the triplon representation are reliable even at 
the mean- field level [3]. Improvements beyond mean- field 
theory can be made by resummation of selected diagrams 
(at lowest order in the density) as typically done in the 
case of a dilute Bose gas [5] . This leads to accuracy within 
10 percent, corresponding roughly to the low quasiparticle 
density n < 0.1. 

It is clear that since for the model Jl} the VBS order is 
spontaneous and is expected to disappear at the QCP, the 
triplons are not necessarily in the dilute limit, especially 
around the QCP. This makes the analysis quite involved 
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Fig. 1. (a.) Lowest triplet (A) and singlet (E s ) energy gaps near 
the QCP. The singlet binding energy e = 2 A — E s . (b.) Triplon 
quasiparticle density n. (c.) Dimer order parameters. 

and the results somewhat dependent on the level of approx- 
imation. Moreover, the destruction of the dimer order at 
the QCP is expected to lead to vanishing of a singlet en- 
ergy scale E s (in addition to the vanishing of the triplon 
gap) . Here we present our main results while the technical 
details will be available elsewhere [6] . 

We have found that the quantum fluctuations in the 
dimer background are strong, which is evident from the 
fact that the four-spin K term in (JTJ favors equally "ver- 
tical" dimer ladders (as in Fig. [2(b)) and horizontal ones, 
obtained by a 90 degree rotation. Since the triplons are 
constructed on a fixed dimer pattern, such a tendency to- 
wards restoration of rotational symmetry is related to the 
tendency of the triplons to form "pairs" with total spin 
S=0, in essence reflecting "plaquettization" . Within our 
approach this is a non-perturbative effect. We denote the 
energy of such a two-triplon singlet excited state by E s (it 
is easy to see that the lowest singlet branch also carries 
zero total momentum), and its wave- function is of the form 
= E a , qH *(q) 4, q 4-q|0). The pairing is most ef- 
fective along the vertical ( "y" ) direction since in this case 
the plaquette interaction acts fully on two dimers, whereas 



pairing involving four triplons is energetically much less 
favorable. The excitation energy w(k) of the one-triplon 
kJO)) s P m S=l state has a minimum at the Neel order- 
ing wave- vector ^af = (ft, ft)', we define the gap as A = 
lo(]laf)- The evolution of A and E s as a function of the 
coupling K/J is shown in Figure 1(a). We find that the 
quantum phase transition to the Neel phase takes place at 
(K/J) c w 2.16, where A — > 0, in good agreement with 
the Monte Carlo result [3] . In addition, the singlet energy 
scale E s also vanishes at the QCP, reflecting the destruc- 
tion of the dimer order. The direct evaluation of the two 
dimer order parameters D x = | (S3 • S4) — (S5 • 84)! and 
D y = I (S3 • S4) — (Si • S4)|, where spins are numbered as 
in Fig. HKb), indeed shows that they exhibit a tendency to 
vanish at the QCP. D x , D y are plotted in Figure 1(c). 

The unconventional merger of singlet and triplet modes 
at the QCP is accompanied however by an increase of the 
triplon quasiparticle density n = ^2 a {tl a tia), as shown in 
Figure 1(b). This signals a breakdown of our quasiparti- 
cle description based on a picture of strongly-interacting, 
but dilute Bose gas of triplons. Within the effective bond 
triplon theory such a breakdown is in a way a necessity if 
the dimer order vanishes at the QCP and the lattice sym- 
metry is restored. In our view the relatively sharp variation 
of all quantities near the QCP is an artifact of our self- 
consistent scheme based on the assumption n <C 1; we thus 
cannot explore reliably the variation of the energy scales as 
(K/J) — > (K/J) c . A consistent classification of diagrams 
beyond the conventionally used "ladder" ones in the dilute 
limit [5] is in itself an unsurmountable task. 

The diverging length scale L s ~ 1/ E s is not the spinon 
deconfinement scale [T] since the state |'3/) s defined above 
does not carry any topological quantum numbers. How- 
ever L s is related to the destruction of the dimer order and 
thus is expected to implicitly reflect deconfinement; L s is 
indeed the most natural scale to introduce within the effec- 
tive triplon theory. The presence of L s in addition to the 
diverging magnetic scale Lm ~ 1/A at the QCP signals a 
tendency towards restoration of lattice symmetry and thus 
supports the notion that the Heisenberg model ([T]) exhibits 
"deconfined criticality" [3] . 
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